SOME NEW SIMPLIFICATIONS FOR GENERAL FORM 
OF ORBIT-ORBIT INTERACTION OPERATOR 



G. Gaigalas 

State Institute of Theoretical Physics and Astronomy, A. Gostauto 12, 

2600 Vilnius, Lithuania 

February 2, 2008 



Abstract 

The general tensorial form of the orbit-orbit interaction operator in the formalism of 
second quantization is presented. Such an expression is needed to calculate both diago- 
nal and off-diagonal matrix elements with respect to configurations, in a uniform way. 
Some special cases are shown for which the orbit-orbit interaction operator reduces to 
simple form. The importance of different contributions to the Breit-Pauli hamiltonian 
is investigated in ground states and ionization potentials for Li isoelectronic sequency 
in the systematic way. 
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1 Introduction 



For light atoms, there are two sorts of corrections to the non-relativistic energies and wave 
functions that frequently are omitted but may need to be included to improve the accuracy: 
the effect of the finite mass of the nucleus and relativistic effects [1, 2]. The lowest order 
corrections for the former, which is particularly important for the very light atoms, can be 
included through a redefinition of the Rydberg constant, Rm = jj+m » f° r converting 
from atomic units to cm -1 , and a mass-polarization correction given by the Hamiltonian 

Kmp = - w S (Pi ■ Pi) > W 

i<J 

where M, m are the mass of the nucleus and the mass of the electron, respectively. Correc- 
tions for the relativistic effects can be obtained by replacing the non-relativistic Hamiltonian 
by the Breit-Pauli Hamiltonian, TCbp, which includes the low-order terms of the Dirac- 
Coulomb-Breit operator, namely terms of the order a 2 (a is the fine structure constant). 
This approach is called as Hartree-Fock-Pauli (HFP) approximation. 

The Breit-Pauli Hamiltonian is often expressed in terms of operators Hi, i = 0, ...,5 
introduced by Bethe and Salpeter [3], but it is also constructive to separate the components 
according to their effect on the spectrum as suggested by Glass and Hibbert [4] , namely 

Hep = 7~Lnr + Hrs + Hfs, (2) 

where TLnr is the ordinary non-relativistic many-electron Hamiltonian. The relativistic 
shift operator Tips commutes with L and S and can be written 

Urs = H-mc + T~(-di + T~(.D2 + T~(.oo + T~(-ssc, (3) 
where TLmc is the mass correction term 

2 N 

Wmc = - T Ep'- ( 4 ) 

The contact interactions describe the one- and two-body Darwin terms Ttoi and 7io2- 
They are: 

Za 2 ir N N 
Hdi = — ^ — 5Z^( r *) and = -na 2 ^25 (rij) . (5) 

i=l i<j 

T~(-ssc is the spin-spin contact term 

87ra 2 N 

T-issc = y~ ^ & ' S ^ S ^ ^ 

i<j 
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and finally Hoo is the orbit-orbit term 
Hoo 



or 
~2 



N 

E 

i<j 



Pj) , ( r ij( r ij -P»)Pj) 



^3 



(7) 



One-particle operators Kmc and Wdi cause relativistic corrections to the total energy. 
Two-particle operators Hd2, Hoo and Hssc define more precisely the energy of each term. 
The Hmc, Hdi, Hd2 and Hssc operators are included into the calculation by adding some 
terms to the radial integrals from non-realitivistic Hamiltonian. The orbit-orbit operator 
needs separate calculations. 



The fine-structure operator Hfs describes interactions between the spin and orbital angular 
momenta of the electrons, and does not commute with L and S but only with the total 
angular momentum J = L + S. So they describe the term splitting (fine structure). The 
fine-structure operator consists of three terms 

Hfs = Hso + Hsoo + Hss- (8) 

The most important of these is the spin-own orbit interaction Hso representing the interac- 
tion of the spin and angular magnetic momentums of an electron in the field of the nucleus. 
The spin-other-orbit TLsoo and spin-spin Hss contributions may be viewed as corrections 
to the nuclear spin-orbit interaction due to the presence of other electrons in the system. 

The two-body terms Hss and Hsoo are complex leading to many radial integrals. The 
complexity of the two-body Hoo operator however, exceeds those of Hfs, increasing the 
computer time required to evaluate an interaction matrix. Thus, it has been customary to 
omit the orbit-orbit effect from energy spectrum calculations. 



The first to treat the Hoo operator in more detail were Dagys et al. [5, 6], and Rudzikas 
et al. [7]. The expressions were obtained for the matrix elements of sp, sd, pp' and dd! 
configurations cases, and for the configurations with a single open shell. Later, the tensorial 
form of this operator and the matrix elements between different configurations were inves- 
tigated by Beck [8], Match and Kern [9], Walker [10], Saxena et al. [11], Dankwort [12], 
Anisimova and Semenov [13], and Anisimova et al. [14]. Among them, Wybourne [15] had 
expressed the matrix elements of this operator inside one open shell through the eigenvalues 
of Casimir operator. The irreducible tensorial form of the orbit-orbit operator was farther 
examined and simplified by Kaniauskas and Rudzikas [16]. 



The general method to calculate the matrix elements of any two-body energy operator 
between arbitrary configurations was proposed in Jucys and Savukynas [17], Eissner et al. 
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[18] (the latter is incorporated into SUPERSTRUCTURE [18]), Glass [19], Glass and Hib- 
bert [4], and was investigated further by Badnell [20]. An interesting recurrent relationship 
for the matrix element of orbit-orbit operator between configurations with one open shell 
was obtained by Kickin and Rudzikas [21]. Exhaustive tables of angular coefficients for this 
operator are presented in Jucys and Savukynas [17]. The matrix elements for configurations 
with /-shells were investigated by Saxena and Malli [22]. 

Investigating the orbit-orbit interaction operator is made a lot easier after rewriting it in 
terms of the products of tensorial operators and then applying the method of Racah. So this 
paper aims to present the orbit-orbit operator in the style Gaigalas et al. [23, 24] and using 
the integral properties find some new simplifications for the operator as well. It enables 
us for evaluation of matrix elements of orbit-orbit interaction operator i) to use full Racah 
algebra, namely allows to calculate both diagonal and off-diagonal matrix elements with 
respect to configurations, in a uniform way; use the unit tensors in three spaces (orbit, spin, 
and quasispin), ii) to take into account new simplifications in general way. 



2 Orbit— orbit interaction 



The tensorial form of orbit-orbit operator (7) in general case is [23]: 

njoo \ ~> ( rr(A:feO,000) . „-(fcfc0,000) . „-(fcfcO,000) . rj(kk0,00)\ 
H =2^ { H ool + ^002 + -^003 + H oo4 ) > 



(9) 



where 



H 



(fcfcO,000) 



ool 



c (*) x c 2 ( fc )l {0) 



„fc-l 

J< 



d d 



2k + 3 rl +2 J dn dr 2 ' 



(10) 



H. 



(fcfcO.OOO) 
oo2 



ia 



k {k + 1) (2k + 1) 



(0) 



k — 2 r. 



2k 



k-l 



k + 1 r\ 



k-l 



k r 



k+l 



— ^T e(n " r2) + 2J—1 Tf e(r2 " ri) + 2TT3 ^p £(ri " r2 



fc + 3 r\ +l \ 1 d 

"2TT^^p e(r2 - ri) n^ (1 + Pl2) ' 



(11) 
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rj(kkO,000) 
U oo3 



-a 2 VW+l 



2k - 1 



fc + 1 



C 



(fc-i) 



x L 



(i) 



C. 



(fc-i) 



x L 



(i) 



(fc) 



(o) _fc-i 



„fc+2 ' 



(12) 



ij(fcfc0,00) 

-"004 



= ^- v / 2fcTi 



C x (fc) x L x 



(i)i (*) 



(1) 



(AO 



(0) 



(A; -2) (fc + 1) r 



fc-2 



fc (fc + 3) 



2k -1 



2fc + 3 



The e(x) in (11) is a Heaviside step-function, 



e(s) 



1; for x > 0, 
0; for x < 0. 



(13) 



(14) 



Here in equations (9) - (13) we imply that a tensorial structure indexed by (fc^fc, o\Oi<j) 
at i? 00 i, -ff o2, #003 and H OQ 4 has rank fci for electron 1, rank k2 for electron 2, and a 
resulting rank k in the I space, and corresponding ranks o\Oiv in the s space. So four 
terms of orbit-orbit operator have the same tensorial structure (kkO 000), summed over k 
in H°° expression (9). It means that orbit-orbit operator is scalar in s space. 



The general expression for any two-particle operator proposed by Gaigalas et al. [24] is 
sutable for evaluation of diagonal and non-diagonal matrix elements in uniform way. It 
allows one to make the most of the advantages of Racah algebra (see Racah [25, 26, 27, 28]). 
So further we will investigate the orbit-orbit interaction in the framework of this formalism. 
This expression has tensorial form: 

iji'j' 



E e(S){4^(n a A a> H)«(«,l) 

a Ki2,ai2,K' 12 ,a' 12 



P,-p 



<5(«,2) 



x D 



(K12CT12) 



E «2<2) (n 7 A 7 ,H) 



(fcfc) 



<5(u,3) 



p,-p 
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[D ilaS) x £>(^ s ) 

/3 7 5 L 



(K12012) 



( K 'l2°12) 



(kk) 



6(u,4)\, (15) 



p-p 



where (i,j\g\i',j') is the two-electron matrix element of operator G^ KlK2k ' <Jia2k \ a\ is the 
electron creation, a}- 
are strictly different. 



electron creation, a}- electron annihilation operators, iji'f, i = niliSimi.m Si and a, j3, 7, 5 



The summation in the (15) runs over the principle and the orbital quantum numbers of 
open shells. The first term represents the case of a two-particle operator acting upon the 
same shell n a \ a , the second term corresponds to operator Qi^^k^^k) ac ti n g U p Qn two 
different shells n a A a , npXp. When operator Q^i^k^i^k) ac ^ g U p 0n three shells the third 
term in (15) must be considered and when it acts upon four - the fourth one. We define 
in this expression the shells n a X a , npXp, n 7 A 7 , n$Xs to be different. In general case the 
number of combinations of iji'j' (distributions) in (15) is infinite. In the work by Gaigalas 
et al. [24] an optimal number of distributions is chosen, which is enough for investigation 
of two-particle operator in general. 



The tensorial part of a two-particle operator in (15) is expressed in terms of operators 
of the type (nX,E), B^ kk \nX,E), C^ kk \nX, £), D^ ls \ E^ kk \nX,E) defined in [24]. 

They denote tensorial products of those creation/annihilation operators that act upon a 
particular electron shell, A = Is, and u is the overall number of shells acted upon by a given 
tensorial product of creation/annihilation operators. Parameter H implies the whole array 
of parameters (and sometimes an internal summation over some of these is implied, as well) 
that connect G with tensorial products of creation/annihilation operators in the expression 
(15) (see [24]). These 6 (H) are all proportional to the submatrix element of a two-particle 
operator g, 

G (H) ~ (niXiUjXj \\g\\ niiXi'rijiXji) . (16) 

So to obtain the general expression of orbit-orbit operator, analogous to (15), the two- 
electron submatrix elements (16) must be defined. In the following section we present these 
explicit expressions for this operator. 



3 Submatrix elements for the orbit— orbit operator 



The sum of submatrix elements of three terms H^ ' 000 ^ , °' 000 - ) and H^ ' 000 ^ is equal 
to (see Badnell [20]): 



( 



rnXirijXj 



H. 



(fcfc0,000) 
ool 



(fefc0,000) 
oo2 



. rr(kk0,000) 
+ H ooi 



Tli' Xi/Tlji Xj/^j 
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= -2[k\ 1 / 2 (li \c {k) \ k>) (lj \\c {k) \\ If) (1 -S(k,0))Z k (niknjlj^idi'nflf) , (17) 
where we have used the conventional shorthand notation [k, ...] = (2k + 1) • ... and where 
Zfc (riiliTijlj ,ni'li>njtlj'^ 

= 2k (k + 1) (r k+1 (mlinjlj,ni>li>njiljf) — T k ~ 1 (riiliTijlj, ni'li'Tij/lj')) 

+ (^ (li + 1) — A; (k + 1) — 1^ (li> + 1)) (ll k+1 (riiliTijlj, rii'li'Tijilj') — U k ~ 1 (riiliTijlj, nidiirijilji)) 
+ (lj (lj + 1) — k (k + 1) — Iji (lj/ + 1)) (ll k+1 (rijljriili, rij'ljirii'li') — U k ~ 1 (rijljnili,njdjiriidii)) 

+ i (k (h + l)-k(k + l)- k, (k, + 1)) (lj (lj + 1) - k (k + 1) - /jv + 1)) 

x - ^ _ yN k ~ 2 (riilinjlj,ni,li,rij,lj,) + N k ~ 2 (rijljUik, nj>lj>ni,li>)j 

k + 3 / \ 
~ (fc+T) (2fc + 3) V fe ( n i l i n j l ji n i' l i' n j' l j') + ^ (njljnili,rij>lj,nidi>)) . (18) 

The radial integrals are defined as 

T fc (riiliTijlj, rii'liirij'lj' 



a 



4(2ife + l) 

i i «C.)^M^t(^ + -)p.-m(^ + -)p,M^,(i 9 ) 



?7 fc (riiliTijlj, ni'li'Tijilji) 



a 



4(2fc + 1) 



/•oo /-oo / ».£.' r k ^ 

/ / ^(n)^-^) (fc-l) nE ^€(r 1 -r 2 )-(fc + 2)-i Tr €(r 2 -r 1 ) 



_/V fe (riiliTijlj , rii'liirij'lj') 



Of roc poo 

= t /„ Jo p * (n) Pj (r2) ^p e(n _ r2)jFV (ri) Pj " (r2) dridr2 ' 

The integrals iV fe (riiliTijlj, Tii'li'Tijilji), T k (riiliTijlj, ni'li'Tijilji) and 
f/ fc (riiliTijlj, ni'li'Tijdj') have the following symmetry properties [12]: 

N k (riiliTijlj , rii'liirij'lj') 

= N k (rii'li'Tijil ji , riiliTijlj) = N k (nidiiTijlj,riiliTijdji) 

= N k (riiliTljilji jTli'li'Tljlj) , 



(20) 



(21) 



(22) 



T k (nilinjlj,ni'liinjilji) = T k (njljndi,njdjmidii) , 



(23) 



U k (nilinjlj,riiiliinjilji) = U k (riidimjlj ,nilinjdj>) 



(24) 



As is seen from the above expressions, the symmetry of the T k and U k integrals is much 
more restricted as compared to the N k integral. There are some useful relations between 
these type of integrals, namely [12]: 

T k (nilinjlj,ni'li'njilji) +T k (riidi'Tijlj ,nilirijilji) = U k (riilirijlj ,nidiirijdji) , (25) 

U k (riilirijlj ,ni/li/rij'ljr) + U k (nidiirijdji ,mlinjlj) 

(jfe — l)(jfc + 2) r i 
= 2k + I {-^ fe_1 (nilinjlj,ni'li>njdj>) + N k ~ 1 (njljnili,nj>lj>mdi') | 

+A (nilinjlj,nidiirijdji) , (26) 

where 

q/2 roo 

A (nilinjlj,riidi>njdj>) = — / Ri (r) Rj (r) R^ (r) Rj> (r) r dr. (27) 

•1 ./() 

We will use these relations in section 4 for getting the simplified expressions for submatrix 
elements of orbit-orbit operator. 

As is seen from the expressions (17) and (18), the matrix elements of H^°' 00 \ H^ ' 00 ^ 
and H^ ' 00 ^ have the same angular dependence as the electrostatic (Coulomb) electron 
interaction operator }{ Coulomb _ s Q it is most convenient to evaluate these three terms 
simultaneously with the electrostatic electron interaction operator [18] which itself contains 
the same tensorial structure 



/ j C oulomb 
k 



and its submatrix element is 



(niXirijXj 



H. 



(fcfcO,000) 
Coulomb 



Hi' Xi'Tlj' Xj'^j 



= 2[fc]V2 (k\\c^\\i tl ) (ijlc^pr) 

n u . • i . c Tj-(kkOfiOO) rj-(/cfcO,000) ^(fcfcO.OOO) , rr(fcfcO,000) ■ 

So submatrix element for H^^^J , H y ool H y oo2 and H y oo4 is 

(n^A, + + ^S°'° 00) + ^f' 000) | n, WA,,) 

x {Rk (nilinjlj,nidi'njdji) — (1 — 5 (k,0)) (rtilirijlj , riidi'iijdj')} . 



(28) 



(29) 



(30) 
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It is more convenient the remaining term H^ o3 ' to calculate separately. Its matrix 
element is 



The use of the approach presented in [24] presumes that both the tensorial structure of 
the operator under consideration and the submatrix elements (niXirijXj \\g\\ ni'Xiirij/Xj/) 
are known. The formulae (17) or (30) and (31) are the expressions we need. We may 
readily obtain the value of a matrix element of this operator for any number of open shells 
in bra and ket functions, by choosing tensorial structure from (9), using their submatrix 
elements in an expression of the type (15), defining bra and ket functions, and performing 
spin-angular integrations according to [24]. 

4 Some simplification for submatrix elements 

In this section we will discuss some special cases of distributions iji'f from Gaigalas et 
al. [24] for the orbit-orbit interaction operator. 

Let us at first consider the distribution iji'f = a[3af3. Using the (18) we express the 
coefficient Z k as 

Z k (n a l a nplp, n a l a nplp) = Z' k (n a l a nplp, n a l a nplp) + Z% (n a l a nplp, n a l a nplp) , (32) 





x (k> -h + k + 1) (lj + lj> + k + 2) x (lj + lj, - k) (lj - If + k + l) 



x (Zj .,_ ii + fc + i))i/2 (y|c(* +1 )||^) (ijWC^Wlj.) 




(31) 



where 



Z' k {n a l a nplp,n a l a nplp) 

= k(k+l) 2T k+1 (nJanplp, n a l a nplp) 

-U k+l {n a l a nplp,n a l a nplp) - U k+1 {nplpn a l a ,nplpn a l a ) 

~ )u k ^\]oul*\ ( Nk ( n «^ V/3> nj a nplp) + N k (nplpr 





(33) 



Z'l (n a l a nplp, n a l a nplp) 
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k(k + l) 2T k 1 {n a l a nplp,n a l a nplp) 



U k 1 {n a l a ni3li3,n a l a nplp) - U k 1 (nplpn a l a ,nplpn a l a ) 



k-i 



fc(fc+l)(fc-2) 
k(2k- 1) 



(34) 



x (iV fc 2 {nj a nplp,n a l a nplp) + N k 2 (nplpn a l a ,nplpnj a f) 

Let us start to evaluate the expression (33). We can rewrite the T k+1 (n a l a nplj3,n a l a rii3li3) 
using the (25) as 

2T k+1 (n a l a nplp,n a l a nplp) 

= T k+l (nj a npl[3, nJanplp) + T k+1 (n a l a nplp, n a l a nplp) 

= U k+1 (nJanplp^JaUplp) . 

With the help of equation (26) we are rewriting the U k+l (nplf3n a l a ,nplpn a l a ) as 

U k+l {nplpn a l a ,nplpn a l a ) 

= \ U k+1 (nplpnda, nplpn a l a ) + U k+1 {nplpn a l a , nplpn a l a ) 
k(k + 3) 



(35) 



N k (nplpn a l a ,nplpn a l a ) + N k (nplpn a l a ,nplpn a l a ) 



2k + 3 

+A {nplpn a l a ,nplpn a l a ) . 

So inserting equations (35) and (36) in the (33) we have: 

, k (k + 1) 

Z' k (n a l a nplp,n a l a nplp) = A{n(3lpn a l a ,nplpn a l a ) 

After similar rearrangements of the expression (34) we have: 
// fe (A; + 1) 



So finally 



Z k {n a l a nplp, n a l a nplp) = 



or 



H. 



(kkOfiOO) 
ool 



+ H 



(A:fcO,000) 
oo2 



(kkOfiOO) 



ooA 



n a \ a 



0. 



(36) 



(37) 



(38) 



(39) 



(40) 



It means that for distributions iji'f = a(3a(3 we do not need to calculate matrix elements 
of the terms #(^°> 000 ) ^ j£(kko,ooo) &n( j ^(fcfco.ooo) ^ a ^ j n & g j m jj ar wa y ^ j g p OSS ible to 

prove that 

Zk (nj a nj a ,nj a nj a ) = Z k {nplpn a l a ,nplpn a l a ) = Z k {nplpn a l a ,n a l a n a l a ) 

= Z k {nd a nplp,n a l a n a l a ) = Z k {nplpnplp,n a l a nplp) = Z k (nplpnplp,n a l a nplp) 
= Z k (nplpn^lj^alan^lj) = Z k {nplpn 1 l 1 ,n a l a n 1 l 1 ) = 0. (41) 
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So for the distributions aaaa, a(3a(3, (3a(3a, f3aaa, a(3aa, (3(3(3a, (3f3af3, (3^a^, 7/^70 we 
do not need to calculate matrix elements of f^fco.ooo) j jj-{kko,ooo) ^(fefco,ooo) ^ ermg ^ ^ OG 

In these cases the orbit-orbit interaction operator contains the term H^ ' 000 ^ only. The ma- 
trix element of this term has the radial integral of only one type, i.e. iV^ -1 (rijljUili, njiljifii'li') 

It is very well known in the literature [17] that the matrix elements of the orbit-orbit 
operator Hoo, ( ns2 1 S\\7ioo\\ n s 2 1 S) and (ns n's 1 S\\Hoo\\ n s n ' s 1 S), are zeroes. It is 
possible to generalize these statements using the results of present paper. We see that for 
direct part of any diagonal matrix elements or the off-diagonal matrix of the type 

{...nl N ...n'l' N ' ... LS\\Hoo\\-nl N±1 .■■n'l' N '^ 1 ... L'S') we need to calculate the matrix ele- 
ment of H^ ' 000 ^ operator only. Using the fact that ^0| IC^ 1 ^ 1 10^ = 0, we strightforwardly 
from (31) find values of these matrix elements in the case I, I' = 0. These values and values 
of exchange part of diagonal matrix elements are equal to zero in this case. This is valid 
for matrix elements between functions with any number of open electron shells. 

Remaining 33 distributions from Table 1 of Gaigalas et al. [24] have all terms H^ k0 ' 000 \ 
^oo2°' 00 °' ) 1 -£foo3°' 00 °' ) an d -ffoo4°' 00 °' ) • For calculation of matrix elements of these distributions 
we need to find the values of T fc±1 (n^n^, nj'/j'rej//j'), U k±1 (njljnili,njdj'ni/li/), 

N k_1 (njljnili,nj'lj'ni'li'), N k ~ 2 (njljnili,njdjinidii) and N k (njljnili,njdj<nidi>) integrals 
(see (18) and (31)). 



5 The effect of the orbit— orbit interaction on ground states 
in light atoms 

Taking into account the relativistic corrections in the Breit-Pauli approximation in the 
configuration interaction method (CI), it is important to know the matrix elements of op- 
erators considered. As was shown in the Section 4 the matrix elements of the orbit-orbit 
operator Hoo, ( ns<2 1 S\\Hoo\\ns 2 1 S) and (ns n's 1 S'| |7^oo| l 7 ^* n's 1 S), are zeroes. There- 
fore, in this approximation for the configuration Is 2 1 S the corrections due to orbit-orbit 
operator appear through the diagonal matrix elements of the types (np 2 1 S\ \Ttoo\ \ n P 2 1 S), 
(nd 2 1 S\\TLoo\\nd 2 1 S), etc., and through the off-diagonal matrix elements. In investigat- 
ing the level ls 2 ns 2 S, the matrix elements (ns 2 n's 2 S\ \Tloo\ \ns 2 n's 2 S), of the orbit-orbit 
operator are equal to zero, too. So the orbit-orbit operator corrections appear through 
the remaining matrix elements where the resulting terms of bra and ket functions coincide. 
Therefore it is plausible that these corrections are unimportant to the absolute values of the 
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level Is 2 1 S'. Of course, one has to investigate into their exact contribution, as compared 
to other relativistic corrections, in aiming at high accuracy of the results. 

5.1 The MCHF method with Breit— Pauli and mass— polarization correc- 
tions 

The computational method for including nuclear and relativistic effects has been described 
in detail elsewhere [2]. Briefly, the wave function ^{jLS) for an atomic state labelled by 
the configuration 7, and term LS is approximated by a linear combination of configuration 
state functions (CSFs), 

M 

*{ n LS) = Y,Ci*{nLS). (42) 

i=i 

Each <&(~fiLS) is constructed from one-electron spin-orbitals for the configuration 7^ and is 
of the same LS symmetry as the atomic state function. In the MCHF method, the radial 
functions used to construct the CSFs and the expansion coefficients c% are determined 
variationally so as to leave the non-relativistic energy stationary with respect to variations 
in the radial functions and the expansion coefficients [2]. Once radial functions have been 
determined, they may be used as a basis for a wave function expansion including additional 
effects. In particular, when relativistic corrections are included, 

M LS 

*( 7 LSJ) = c^lsH^LSJ). (43) 

LS i=l 

where ^(jLSJ) is a CSF in which the total angular momentum L and the total spin S have 
been coupled to a resultant J. The expansion coefficients are obtained from a configuration 
interaction calculation, where the interaction matrix is evaluated with respect to the Breit- 
Pauli Hamiltonian, or some subset of operators. New, efficient programs based on the 
combination of second quantization in coupled tensorial form, and a generalized graphical 
technique [24] were used for evaluating the Breit-Pauli operators. The expressions for 
orbit-orbit interaction operator are taken from present work. 

5.2 Ground states and ionization potentials for Li-like atoms and ions 

The ground states of Li isoelectronic sequency was found using method described above. The 
calculations reported here are strictly ab initio: no ^-extrapolation or basis extrapolation 
has been applied. The configuration states included in the expansions of different terms were 
obtained by including all possible CSFs of a given LS symmetry that could be constructed 
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Table 1: Comparison of contributions to the Breit-Pauli energies (in au) between MCHF 
(present work) and the full core plus correlation results of Chung [29] (the second line). 
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Table 1: (continued) 
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from orbitals with n < 10, I < 7. The largest expansion for this rule-based scheme was for 
ls 2 2s 2 S, where the interaction matrix size was 7 496. 

In Table 1, we compare contributions to the Breit-Pauli energies (without the mass- 
polarization correction) of ls 2 2s 2 S and Is 2 x 5 with those reported by Chung [29]. We see 
that the non-relativistic energies are in close agreement though, the present results, without 
any extrapolations are slightly lower than those of Chung. The non-relativistic ionization 
potential (IP) is in close agreement. The relativistic shift effect is in surprisingly large 
disagreement for B III - F VII though, again, the contribution to IP agrees more closely. 
Finally, the much smaller orbit-orbit effect is in good agreement, with the difference again, 
agreeing to more decimal places than the individual energies. Some differences are expected 
since, in the present work, effects are included in the interaction matrix, prior to matrix 
diagonalization whereas in the full-core plus correlation method employed by Chung, these 
result are computed as a small perturbative correction from the non-relativistic wave func- 
tion. However, in the present methodology, there also may be basis effects, in that the 
orbitals used in the expansion are optimized for the non-relativistic Hamiltonian and are 
incomplete with respect to the Breit-Pauli Hamiltonian. This is particularly true for the 
relativistic shift without orbit-orbit interaction where the one-electron Darwin term depends 
only on the value of the s-orbitals at the nucleus. The effect of orbit-orbit interaction on 
the energies is small, and the contribution to the ionization potential is in good agreement 
between the two theories. Finally, to gain some perspective on the magnitude of corrections, 
we include the QED correction of —0.0001344 au. to the ionization potential of O VI as 
reported by Chung [30]. However, it should be noted that the orbit-orbit correction to the 
2p ionization potential is -0.0003328 au. [31] so the relative importance depends on the 
state. For a correct spectrum, it appears that both should be included, at least for the 
lower levels. 

In Table 2 and Table 3, we report the contributions to the Breit-Pauli ground energies of 
F VII and F VIII and ionization potential using seven different expansions (first column 
of Table 2 and Table 3). The notation nl < 3d, for example, means that the expansion 
was obtained by including all possible CSFs of a given LS symmetry that could be con- 
structed from orbitals with n < 3, I < 2. It contains 27 configurations for ls 2 2s 2 S and 10 
configurations for Is 2 1 S (see Table 3). 

The Table 2 indicates that in case the CSF number is being increased, non-relativistic 
energy E nr lowers. The energy is also being lowered by Ers-oo and E OQ corrections. Their 
absolute values increase insignificantly at the increase of CSF number. Having compared 
those results with Chung results, we notice that the values of Ers-oo and E OQ indicated in 
the article tally best with the values of Chung, when CSF are generated from orbital nl < 3d 
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Table 2: Comparison of E nr , Ers- 00 , E oq energies for F VII and F VIII ground states in 
seven different expansions. 
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Table 3: Comparison of ionization potential (in au) for the ls 2 2s 2 S states of F VII in 
different expensions with different corrections. IP nr - nonrelativistic case, IPrs-oo ~ only 
relativistic shift without orbit-orbit term, IP 00 ~ only orbit-orbit term. 
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or nl < 4/. Whereas non-relativistic energy values, got in the article, tally best, when CSF 
basis includes orbitals nl < 91. The same discussions valid for ionization potentials in 
different approximations (see Table 3), too. 



6 Conclusions 



The general irreducible tensorial form of the orbit-orbit interaction operator in the for- 
malism of second quantization is presented (expressions (15), (17) (31) ). It contains four 
different terms. Each is associated with different set of radial integrals. In the present work 
we have succeeded in obtaining simpler expressions having only the term H^°' 000 \ for 
some special distributions of electrons in the configuration. As we see from the paper the 
formalism for evaluation of matrix elements developed by Gaigalas et al. [24] allow us to use 
these simplifications for practical applications in general way. This facilitates practical cal- 
culations of matrix elements without restraining the generality, and is one more advantage 
of the approach used. The properties (39) and (41) are useful for testing the calculation of 
matrix elements and for evaluation of the accuracy of radial integrals T k , U k and N k , too. 

The results from section 5 shows that orbit-orbit operator and QED corrections may be 
of comparable size in light elements, though the orbit-orbit operators is of order 0(a 2 ) 
and QED corrections 0(a 3 ). It related with the fact that on the one hand a big number 
diagonal matrix elements of orbit-orbit operator are zeroes, on the other hand off-diagonal 
matrix elements of orbit-orbit operator is of order C(a 4 ). So, need to take into account 
bouth orbit-orbit operator and QED corrections for studing the Li-like sequence and other 
light elements. 
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" Supaprastintos orbit-orbit saveikos operatoriaus formos taikymas lengviems atomams." 

G. Gaigalas 

Santrauka 

Straipsnyje pateikta bendra orbit-orbit saveikos operatoriaus tenzorine forma. Ji susideda 

, , • , • , • ■ Tj(kkO,000) Tr(kkO,000) Tj(kkO,000) ■ Tj(kkO,000) , . ~. . 

is keturiu skirtmgu narm H y ool , H y oo2 , H 00 ^ ' lr H y ooA , turinciu ta pacia 

tenzorin§ struktura, taciau skirtingas radialiasias dalis. Autorius surado tokius atvejus, kai 

darbe nagrinejamas orbit-orbit saveikos operatorius issireiskia per viena nari H^ ' 000 ^ t.y. 

igyja zymiai paprastesn§ israiska nei buvo zinoma iki siol. Straipsnyje pasiulytas budas, 

kaip sias naujas israiskas butu galima naudoti bendrai, t.y nepriklausomai nuo to tarp kokiu 

konfiguraciju ieskomi sio operatoriaus matriciniai elementai. Tai leidzia: i) tiek diagonalius 

tiek nediagonalius (konfiguraciju atzvilgiu) matricinius elementus nagrineti vieningai, ii) 

efektyviau algoritmizuoti orbit-orbit saveikos operatoriaus matriciniu elementu skaiciavima, 

iii) atlikti radialiujy integralu tikslumo ivertinima. 

Remiantis straipsnyje pasiulyto metodo pagrindu atlikti teoriniai Li izoelektronines sekos 
pagrindines btisenos ir ionizacijos potencialo skaiciavimai. Siuo pasirinktu atveju, atsiranda 
nemazas skaicius tokiu matriciniu elementu, kuriuos nagrinejant naudojamos supaprastintos 
orbit-orbit operatoriaus israiskos. Darbe gauti teoriniai rezultatai sutapo su teoriniais kitu 
autoriu rezultatais. Tai itikinamai parodo, jog darbe gautos supaprastintos israiskos bei 
pasiulyta metodika (kaip jas efektyviai isnaudoti) yra teisinga. 

Nagrinejant sudetingesnius atomus bei jonus, orbit-orbit operatorius ne visada turi su- 
paprastinta pavidala. Tuo atveju si operatoriu tikslinga nagrineti kompleksiskai t.y. kur 
imanoma naudotis supaprastintoms, kur ne bendrosiomis israiskomis. Si problema taip pat 
straipsnyje isspr§sta. 
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